WELL-POSEDNESS OF AN INTEGRO-DIFFERENTIAL 
EQUATION WITH POSITIVE TYPE KERNELS MODELING 
"(N ■ FRACTIONAL ORDER VISCOELASTICITY 

o 

CSl ■ FARDIN SAEDPANAH 

Abstract. A hyperbolic type integro-differential equation with two weakly 
singular kernels is considered together with mixed homogeneous Dirichlet and 
non-homogeneous Neumann boundary conditions. Existence and uniqueness 
of the solution is proved by means of Galerkin's method. Regularity estimates 
are proved and the limitations of the regularity are discussed. The approach 
presented here is also used to prove regularity of any order for models with 
smooth kernels, that arise in the theory of linear viscoelasticity, under the 
appropriate assumptions on data. 
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< 

3 . 1- Introduction 



We study the model problem (|2.6|) , which is a hyperbolic type integro-differential 
equation with two weakly singular kernels of Mittag-Leffler type. This problem 
arises as a model for fractional order viscoelasticity. The fractional order viscoelastic 
model, that is, the linear viscoelastic model with fractional order operators in the 
f~^ ■ constitutive equations, is capable of describing the behavior of many viscoelastic 

f^ I materials by using only a few parameters. 

^^ ■ A perfectly elastic material does not exist since in reality: inelasticity is always 

fT^ I present. This inelasticity leads to energy dissipation or damping. Therefore, for 

^^ i a wide class of materials it is not sufficient to use an elastic constitutive model 

Cn I to capture the mechanical behaviour. In order to replace extensive experimental 

tests by numerical simulations there is a need for an accurate material model. 

Therefore viscoelastic constitutive models have frequently been used to simulate the 

k> ■ time dependent behaviour of polymeric materials. The classical linear viscoelastic 

'^ I models that use integer order time derivatives in the constitutive laws, require an 

C^ I excessive number of parameters to accurately predict observed material behaviour. 

Bagley and Torvik [3] used fractional derivatives to construct stress-strain rela- 
tionships for viscoelastic materials. The advantage of this approach is that very 
few empirical parameters are required. When this fractional derivative model of 
viscoelasticity is incorporated directly into the structural equations a time differen- 
tial equation of non-integer order higher than two is obtained. One consequence of 
this is that initial conditions of fractional order higher than one are required. The 
problems with initial conditions of fractional order have been discussed by Enelund 
and Olsson j^, see also references therein. To avoid the difficulties with fractional 
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2 F. SAEDPANAH 

order initial conditions some alternative formulations of the fractional derivative 
viscoelastic model are used in structural modeling. The formulation that we use, is 
based on a convolution integral formulation with a singular kernel of Mittag-LefSer 
type, see [I], [8], and [9]. 

There is an extensive literature regarding well-posedness and numerical treat- 
ment for integro-differential equations, see, e.g., [1], [4], [11], [12], [13], [H], [15], 
[16], [17], and [18]. Existence, uniqueness and regularity of a parabolic type integro- 
differential equation has been studied in [2] by means of Fourier series. One may 
also see [7], where the theory of analytic semigroups is used in terms of inter- 
polation spaces to solve a boundary value problem in linear viscoelasticity. An 
abstract Voltcrra equation, as an abstract model for equations of linear viscoelas- 
ticity, has been studied in [F . In a previous work [llj , well-posedness and regularity 
of the problem ()2.1ip . which is a simplified form (synchronous viscoelasticity) of the 
model problem (j2.6p . was studied in the framework of the semigroup of linear opera- 
tors. The drawback of the framework is that this does not admit no n- homogeneous 
Neumann boundary condition, while in practice mixed homogeneous Dirichlet and 
non- homogeneous Neumann boundary conditions are of special interest. Here we 
investigate existence, uniqueness and regularity of the solution of the problem ()2.6|) 
by means of the Galerkin approximation method. 

There are also models with exponential kernels, smooth kernels, which describe 
polymeric materials, e.g., natural and synthetic rubber. Existence, uniqueness, 
and regularity of the solution of such models can be adapted from, e.g., [5]. The 
drawback with this kind of models is that it requires a large number of exponential 
kernels to describe the behavior of the materials. This is the reason for introducing 
kernels of Mittag-Leffler type or fractional operators. In [2 and 9 it is shown 
that the classical viscoelastic model based on exponential kernels can describe the 
same viscoelastic behaviour as the fractional model if the number of kernels tend 
to infinity. Here we extend the presented approach so that regularity of any order 
of the solution for the models with smooth kernels can be proved. 

In the sequel, in §2 we describe the construction of the model problem (|2.6p and 
we define a weak (generalized) solution. Then in §3, using Galerkin's method, we 
prove existence and uniqueness of the weak solution of the problem. Finally, in §4 
we study regularity of the solution and limitations for higher regularity. We also 
show that higher regularity of any order of the solution of models with smooth 
kernels can be achieved. 

2. Fractional order viscoelasticity 

Let cr.y , Cij and Ui denote, respectively, the usual stress tensor, strain tensor and 
displacement vector. We recall that the linear strain tensor is defined by, 

1 / dui duj ' 
'■^ 2 V dxj dxi 

With the decompositions 

the constitutive equations are formulated as, see [31, 

s,j(i) + T'^'Dt's.jit) = 2Gooe^,(t) + 2GT'^'Dt'e,j{t), 
^^-'^' akk{t) + T^'D^-'akkit) = 3Kooekk{t) + SKr^'D^'ekkit), 
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with initial conditions 

s,,(0+) = 2Ge,j(0+), afefe(0+) = 3i^efefe(0+), 

meaning that the initial response follows Hooke's elastic law. Here G, K are the 
instantaneous (unrelaxed) modulus, and Goo, K^o are the long-time (relaxed) mod- 
ulus. Note that we have two relaxation times, ri,T2 > 0, and fractional orders of 
differentiation, ai, a2 G (0, 1), where the fractional order derivative is defined by 

D1f{t) = DtD;^'-''^f{t) = Dt—^ f\t - s)-"f{s) ds. 

r(i - a) Jq 



The constitutive equations (|2.ip can be solved for a by means of Laplace transfor- 
mation, ^: 

{t - s)e.,j{s)dsy 
2{t~ s)ekk{s)ds 





S^,{t) = 


^2G(e..,{i)-^^j\, 




CTfefe(t) = 


^34.(0-^/-/. 


where, for i - 


= 1,2, 


d „ / / i \°"\ „ ... 



<"-s^-(-W°')> ^-("^E 



„3,r(i + „<,.)• 

and i?Q; is the Mittag-Leffler function of order ai. Then we define parameters 7^, 
and the Lame constants /i. A, 

G — GcxD -fi ^ jVp, 

"g ' ^' " iT 

We also define Pi — jiOi, and the constitutive equations become 



71 = — -77-^, 72 == "' rx"^°° ' A* = G, A = /-i:-|G. 



(T„(i) = \2fieij{t) + Xekk{t)5ij] ~ 2^ /?i(t - s)Uij(s) - ■|efcfc(s)%j ds 

3A-t-2A. /■* 
5 / P2[t- s)ekk{s)d^jds. 

«J JO 

The kernels are weakly singular, for i = 1, 2: 
and we note the properties 

P^it) > 0, 
IIAI1li(R+) = y |3^{t)dt - 7.(£^a.(0) - S„.(oo)) = 7, < 1 . 

The equations of motion are 

pui^tt - cTijj = f'i., in fl, 

(2.4) Mj = 0, ouTd, 

ai-jUj ^ gi, on Tn, 

where u is the displacement vector, p is the (constant) mass density, / and g 
represent, respectively, the volume and surface loads. We let fl C M'', d — 2,3, 
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be a bounded polygonal domain with boundary Fd U Fn = dfl, Fd H Fn = and 
meas(FD) ^ 0. We set 

{Au)i = -(2^eij(u) + Xekk{u)Sij) 



J 



(Aim)., = -2/i(e,y(M) - ■iefcfc(M)(5ij) ■ , 



(2.5) i^^ru^., = -z/l^^ty■v"^ - 3 

{A2U)^ = ^ {ekkiu)S^J)_. , 

and clearly we have A ^ Ai + A2. Now, we write the equations of motion p.4p in 
the strong form, (we denote time derivatives with '•'), 

pu{x,t) + Au{x,t) 

2 ,t 



V I /3,{t-s)A,u{x,s)ds = f{x,t) 
.=1 -^0 



in f7 X (0,T), 



(2.6) M(x,t) = onFDx(0,r), 
tT(u; x.t) ■ n ~ g{x, t) on Fn x (0, T), 
m(x, 0) = u^{x) in J7, 

M(a;,0) = f'^(a;) in 51. 

We define the bilinear form (with the usual summation convention) 

a{u,v) ^ / {2fieij{u)eij{v) + Xeu{u)ejj{v)) dx, VM,weV", 

which is coercive, and in a similar way, corresponding to Ai, the bilinear forms 
ai{u,v), i ~ 1, 2. We introduce the function spaces H = L2{^)'^, i?rN = ^2(rN)'^, 
and V = {v ^ H^{i^Y ■ ^|rD= 0}. We denote the norms in H and Hr^^ by || • || 
and II • IItnj respectively, and we equip V with the inner product o(-, •) and norm 
\\v\\y = a{v, v). It can be shown that A, A\ and A^ are self-adjoint, positive definite 
linear operators (not bounded). We also note that, for i; G T^, 

(2.7) a,[v,v)<\\vfy. 

Now we define a weak solution to be a function u ~ u(x, t) that satisfies 

(2.8) u £ L2iiO,T);V), ii e L2iiO,T): H), il e L2{{0,T);V*), 

2 „t 

(2.9) p{u{t),v) + a{u{t),v) ~y2 / Pi{t- s)ai{u{s),v)ds 

= {f{t),v) + i9{t),v)r^ , ■iveV, a.e. t e (0, T), 

(2.10) u(0) = u°, u{0) = v°. 

Here {g{t), v)r^ — /p g{t)-v dS, and (•, •) denotes the pairing of V* and V. We note 
that ()2.8p implies, by a classical result for Sobolev spaces, that u G C([0, T]; iJ), li G 
C([0,T]; F*) so that the initial conditions (P1I7| make sense for u° G i/, w^ G y*. 

Remark 1. If wc make the simplifying assumption (synchronous viscoelasticity) : 

a ^ ai = a2, r == ri == T2, 9 ^ 9i ^ 02, 

we may define 7 = 71 = 72, so that 13 = Pi = /?2- Then the equations of motion is 

(2.11) pu{x, t) + Au{x, t)- I (3(t - s)Au{x, s) ds = f{x, t) in fl x (0, T), 
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together with the boundary and initial conditions in ()2.6p . 

Well-posedness and regularity of the simplified problem ()2.1ip was studied in 
[11) . in the framework of the semigroup of linear operators. The drawback of 
the framework is that this does not admit non-homogeneous Neumann boundary 
condition. 

3. Existence and uniqueness 

In this section we prove existence and uniqueness of a weak solution of (|2.6p 
using Galerkin's method, in a similar way for hyperbolic PDE's in [5], [TU]. To this 
end, we first introduce the Galerkin approximation of a weak solution of (J2.6I) in 
a classical way, and we obtain a priori estimates for approximate solutions. These 
will be used to construct a weak solution, and then uniqueness will be verified. 

We recall (|2.2[) . (|2.3p and we define the functions 

i^{t) = 7* - / /3^(s) ds^ /3,{s) ds = J^EaM, « = 1, 2, 

Jo Jt 

and it is easy to see that 

(3.1) AC.W = -/3.(i) < 0, 6(0) =7», lime^(i)=0, 0<Ut)<l^■ 
Besides, ^i are completely monotonic functions, that is, 

i-iyDU,it)>0, ie(0,oo),jeN, 

since the Mittag-Leffler functions £"„ . , a^ G [0, 1] are completely monotonic. Conse- 
quently, an important property of 6, i = 1, 2, is that, they are positive type kernels, 
that is, they are continuous and, for any T > 0, satisfy 

(3.2) / f C^{t^s)(f>{t)(l){s)dsdt>0, V</>eC([0,r]). 



3.1. Galerkin approximations. Let {{^j,'Pj)}7!Li be the eigenpairs of the weak 
eigenvalue problem 

(3.3) a{ip,v)^X{ip,v), \fveV. 

It is known that {(pj}°Zi can be chosen to be an ON-basis in H and an orthogonal 
basis for V. 

Now, for a fixed positive integer to G N, we seek a function of the form 

771 

(3.4) u^{t)=J2dj{t)^j 
to satisfy 



(3.5) 



,-1 "'0 



7 = 1 



- (/(t), ifik) + ig{t),ipk)r^, fc = 1, . . . , TO, t e (0, T), 
with initial conditions 

m m 

(3.6) U^(0) =^(w",(pj)<^j, Um(0) =^(w'',(pj)<^j. 
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Theorem 1. For each m E N, there exists a unique function Um of the form ()3.4p 
satisfying dSU-dUl). Moreover, if u° e V,v° e H, g e Wl{{0,T); Hr^), f £ 
L2{{0,T); H), there is a constant C — C{Q, ^1,^2, PtT) such that, 

kni|Uoo((0,T):y) + ||Wto||l^((0,T);H) + ||"m ||l2((0,T);V*) 
^^■^^ < C{\\u'\\v + \\v"\\ + \\9\\wlao,T),Hr,) + II/IIl.((O.T);//) } ■ 

Proof. Using p.4p and the fact that {'fij}°^i is an ON-basis for H and a solution 
of the eigenvalue problem (|3.3p . we obtain from (|3.5I) that, 

m 2 

pdk{t) + \kdk{t) -'^'^ai{ipj,ipk){Pi *dj){t) 

= fk{t) + gk{t), fc = l, ...,TO, te (0,r), 

where * denotes the convolution, and fk{t) = {f{t),ipk), gk{t) = (5(^^)1 Vfe)rN- This 
is a linear system of second order ODE's with initial conditions 

(3.9) dfc(0) = (u°,v'fc), dfc(0) = (i;°,(^fe), fc = l,...,m. 

The Laplace transform can be used, for example, to find the unique solution of the 
system. 

Taking the Laplace transform of p.Sp we get (we use an over-hat for the Laplace 
transform) , 

m 2 

[s^ p -V \k)dk{s) ~ y^y^a,(y)j,<y9fc)^,(s)rfj(s) 

^^■^^^ = Ms) + gk{s) + spdkiO) + pdk{0), 

k — 1, . . . ,m, s > 



rin{ri,T2}' 

Now we show that there is a unique solution for this system, that is, the inverse 
Laplace transform is computable. This proves that the linear system (|3.8p with 
initial conditions (13. 9p has a unique solution, that will prove the first part of the 
theorem. 

Indeed, the Laplace transform of the Mittag-Lefflcr function is 

dE^iaf^)) ^ ^^ , s>|a|i/". 
s"' — a 

Hence for the kernels Pi, i — 1,2, defined in (12. 2p we have 

/3.(s) = Cmt)) = -j,sC{Ea,X-rr''H"')) +j,Ea,M 



^' <1, s>Tr\ 



{nsY^ + 1 



Since Pi are decreasing functions, at most at one point, say s = s > 1/ min{ri, T2}, 
we might have /3i(s) = /?2(s)- Therefore without loss of generality we can assume 
Pi > /32 for s > s. Then equation p.lOp . using a(-, •) = ai(-, •) + a2(-, •)) can be 
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written in the form 

(Ci(s) + (/3i(s) - ^2(s))A2)D(s) = F(s), if s > S, 

(C2(s) + (/32(s)-/3i(s))Ai)D(s) = F(s), if ^^^ r<s<S, 

min{Ti,T2} 

where 

D(s) = (4(s))r=i, F(s) - (A.(s) + .gfe(s) + spdkiO) + pdkiO))Zv 

A, = {a,{<fij,ipk))J,^^-^, C,(s) =diag{(sV + Afc - Afc/3j(s))"^J, j = l,2. 

Obviously C^ and A^, i = 1, 2 are symmetric positive definite matrices, and there- 
fore the coefficient matrices of the finear systems above are symmetric positive 
definite. That is, they are uniquely solvable and this proves the first part of the 
theorem. 

Now we prove the a priori estimate p.7p . Since f3i{t — s) ~ DsCi{t ~ s) and 
Ci(0) = 7i J by p.ip . we can write p.Sp . after partial integration in time, as 

2 

p{u,nit),ipk) + a{ujn{t),(pk) ~y^^liatiu„i{t),ipk) 

2 „t 

+ y2 £.i{t - s)ai{u,n{s),(pk) ds 

2 

-^^i{t)ai(u^{0),ipk), fc = 1, ..., m, i e (0,r). 



Then multiplying by dk(t), summing over k — 1, . . . ,m^ and integrating with respect 
to t, we have 

2 ^t pr 

p||um(i)r+(l -7)ll'"m(t)||y + 2^ / / £,i{r ~ s)ai{um{s),Um{r))dsdr 

<p||u„(0)f + (l-7)||w,„(0)||2. 

+ 2 / (/(r), um{r)) dr + 2 [ {g{r), u„(r))r„ dr 
Jo Jo 

2 ^t 

-2^ / (i{r)ai{u„i{0),Um{r)) dr, 



where 7 = max{7i, 72} and 7 = min{7i, 72}. We note that < 7,7 < 1. Since 
^i, i = 1, 2 are positive type kernels, recalling p.2p . the third term of the left hand 
side is non-negative. Then integration by parts in the last two terms at the right 
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side yields 

p\\Ura{t)f + {^ - i)\WUt)fv 

< p||u™(0)f + (1 - 7)ll"m(0)||y + 2 / {f{r),ura{r))dr 

Jo 

- 2 / {g{r), u™(r))r« dr + 2{g{t),u.^{t))r^ - 2(5(0), M„(0))r« 
Jo 
2 .t 
-2^ / /3i(r)ai(M„i(0),M„i(r))(ir 

2 2 

-2^e^(i)a^(um(0),^in^(i)) + 2X!^^(0)a»("™(0)'"™(0))• 

i=l 4=1 

This, using the Cauchy-Schwarz inequahty, the trace theorem, ||/3i||ij(R+) = 7i, 
S,^{t) < 6(0) = lu and dUT]), imphes 



p||u™(t)f + (l-7)||ti„(t)||2, 

<p||u„(0)f + (l-7)||w„(0)||2. 



+ 2/C1 max ||u„,(r)f + Ci( / \\f{r)\\dr) 

0<r<t \Jq / 

+ 2CTraco/C2 max ||Um(r)||y +2CTracoC2( / \\g{r)\\H^n] 
0<r<t \ /n / 




2 I nri r< II„/-+M|2 



+ 2CTracc/C3||u,„(0lly + 2CTraccC3 || <?(f ) 1||,^ 
+ 2CTracc/C4||"m(0)|i^ + 2C'TraceC4 ||5(0) || 1,^ 



2/^5 ( ^ 7^) l|um(0) 11^ + 2C5 ( ^ 7i) max ||u„(r 

2 2 2 

2/C6(^7.)lkm(0)fv- + 2C6(5]7.)ll«mWlly + 2(5]7.)ll"m(0)| 



1=1 



Hence, considering the facts that Cxraco = C{^), ||um(0)|| < ||w°||, and ||um(0)||y < 
||M°||y, for some constant C = C(r2, 71, 72, /3,r), we have 



l"™llioo((0,T);H) + ll"nillLoc((0,T);V) 



|2 I ||A||2 



<C{||«°f + ||z.°f^ + ||g||i^((o,^)„ 



+ II/IIli((0,T);H) + ll5llLi((0,T);ifrN)}- 

This, and the facts that \\g\\L^((o.T)-HT^) < C\\9\\wiao,T)-Hr^)' by Sobolev inequal- 
ity, and ||/|Ui((o,T);ff) < C'll/l!L2((o,T);if), imply 

ll"™llioo((0,T);i/) + lkm||L<,„((0,T);l^) 
^^■^^^ < C{\\VY + \\uYv + \\9rwlao.Ty,Hr^) + Il/IIL((0,T);«)}- 

Now we need to find a bound for ?i.„i- For any fixed v (£ V with ||u||v' < 1, we 
write V = v^ + w^, where v^ G spa.n{Lpj}^i, v"^ G span({(^j}^]^)^. We note that 
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■y^llv < 1- Then from p.Sp we obtain, 

p{u,n{t),v) = p{um{t), v^) = {f{t), u^) + {g{t), v'^)r^ - a{um{t), w^) 

2 „t 






that, using the Cauchy-Schwarz inequahty, the trace theorem, and (J2.7I) . imphes 

-, 2 

(u™(i),w)| < -f||/(i)||+CTraco||3W|| + ||w„,(i)||y+ max ||u™(s)||y Vt. 

^ i=l 

This, using (J3.12I) . in a standard way imphes 

ll^m|lL2((0,T);V*) ^ '^'{\\f\\L2{iO,T)-H) + \\9\\L2i{0,T)-Hr^) 

+ \\-Y + \\uYv + \\9\\lvimTy,HrJ- 
Therefore, for some constant C — C(ri, 71, 72, p,r), 

ll"m|lL((0,T);y*) < C{\\vY + \\uYv + hWlvimTy.Hr^) + ll/llL((0,T);ff ) } • 

This and (13.12^ imply the estimate p.7p , and the proof is complete. D 

Remark 2. For the simplified model problem (|2.1ip in Remark 1, since /3i = (32, 
the linear system p. lip is reduced to 

CD = r, 

where C = Ci = C2- This system has a unique solution, since C is symmetric 
positive definite. Hence there exists a unique function Um of the form (j3.4p satisfying 

Jo 

- U{t).Vk) + {git), ipk)r^ , fc = 1, . . . , m, i G (0, T), 

with initial conditions p.6p . Moreover, the a priori estimate p.7p still holds with 

3.2. Existence and uniqueness of the weak solution. Now, we use Theorem 
1 to prove existence and uniqueness of the weak solution of (12.61) , that is a solution 

of d^ - firm . 

Theorem 2. If u° G V, v° € H, g e Wl{{0,T);Hr^), f e i2((0,T); iJ), there 
exists a unique weak solution of (j2.6p . 



Proof. 1. We note that the estimate (|3.7p does not depend on m, so we have 

lkm||Lcx,((0,T);y) + ll"m|U,o((O.T);H) + W^mW L2{{0,T)-V) 

<X = if(r!,7i,72,T,uO,t;0,/,5). 
That is, 

{um}T is bounded in LooiiO,T);V) C L2((0,T); 1^), 
(3.13) {um}T is bounded in L^{{0,T)-H) C L2{{0,T);H), 

{ilraJT is bounded in L2{{0,T);V*). 

2. First we prove existence. From p.l3p and a classical result in functional anal- 
ysis, we conclude that the sequencess {um}m=i^ {'^m}m=i^ {^m}m=i ^^^ weakly 
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precompact. That is, there are subsequences of {um}m=i: {'^'"}m=i; {^m}m=iJ 
such that 

ui^u m L2iiO,Ty,V), 

(3.14) ui-^it in i2((0, T);iJ), 

ili^u in i2((0,T);V^*), 

where the index / is a replacement of the label of the subsequences and '^' denotes 
weak convergence. Consequently, (|2.8p holds true and we need to verify (|2.9p and 
p.lO|) . To show (|2^ we fix a positive integer A^ and we choose v e C([0, T]; F) of 
the form 

N 

(3.15) «(i)-^/^,(i)¥',. 



Then we take I > N and by (|3.5p we have 

cT , 2 „t 



fp(iti,v) +a(u;,v) -^ / fii{t-s)ai{ui{s),v)ds)dt 

= / ((/,«) + (g,i;)r.)rfi. 
Jo 

This, by p.l4p . implies in the limit 

/ (p{u,v)+a{u,v)-22 I3i{t - s)ai{u{s),v)ds] dt 
(3.17) -^^ -1-^° 

{{f,v) + {g,v)r^)dt. 



Since functions of the form p.lSp are dense in L2{{0,T);V), this equality then 
holds for all functions v £ L2((0,T); V), and further it implies p.9p . 

Now, we need to show that u satisfies the initial conditions (|2.10p . Let v £ 
C^([0, T]; F) be any function with v{T) = v{T) — 0. Then by partial integration in 
p.l6p we have 



T 





(p{ui,v)+a{ui,v)-^^ / /3i{t - s)ai{ui{s),v)ds] dt 

((/, i;) + (5, «)r„) di - p{ui{0), i>(0)) + p(uz(0), «(0)), 
so that, recalling p.l4p and p.6p . in the limit we conclude, 

ip{u,v)+a[u,v)~y / l3i{t — s)ai{u{s),v) dsj dt 



i=i ■ 
((/, ^') + (5, «)r.) dt - p(u", {,(0)) + piv°, z;(0)). 
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On the other hand integration by parts in p.l7p gives, 



i=l 
T 



(p(u,v)+a{u,v)-22 / l3i{t - s)ai{u{s),v) dsj 



dt 



((/, v) + (g, v)r^) dt - p{u{0),v{0)) + p{v{0),v{0)). 

Conipairing the last two identities we conchide (j2.10p . since v{0), v{0) are arbitrary. 
Hence m is a weak solution of (|2.6p . 

3. It remains to prove uniqueness. To this end, we show that u = is the 
solution of (E^-dOni) ioT u° = v° = f = g = 0. Let us fix r G [0, T] and define 



, , u(u}) duj < t < r 
^0 r <t<T. 

We note that 

(3.18) v{t) e V, v{r) = 0, v{t) = -u{t). 

Then inserting v in (|2.9p and integrating with respect to t, we have 

(3.19) / {p{u,v)+a{u,v))dt~^ l3,{t - s)a,{u{s),v{t)) dsdt = 0. 
Jo i^i Jo Jo 

For the second term, recalling —l3i{t) — Dt^i{t), i — 1,2 from p.ip . we obtain 
l3,{t- s)a^{u{s),v{t))dsdt= / / Dt^,{t - s)a^{u{s),v{t)) dtds 

Jo Js 

^i{r - s)ai{u{s),v{r))ds 


S,i{0)ai{u{s),v{s))ds 



^i{t - s)ai{u{s),v{t))dtds 

= ~Ji / ai{u{s),v{s))ds 
Jo 

+ / / ^i{t ~ s)ai{u{s),u(t))dsdt, 
Jo Jo 

where we changed the order of integrals and we used integration by parts, ^i(O) = ji 
from (13. ip . and v{r) = from p. 181) . Therefore integration by parts in the first 
term of (|3.19l) yields 



[— p{u,v) + a{u,v)) dt — y 7i / ai{u,v)dt 

«=i -^0 

^ pr pt 

+ ^ / ^,{t-s)a^{u{s),u{t))dsdt^O. 



i=l 



_i ^0 JO 



12 F. SAEDPANAH 

This, using p.lSp . implies 

2 

p||^.(r)|p_p||^,(0)f-||^;(r)f^ + |l«(0)f^+^7.(a.Mr),^;(r))-a,(^;(0),^;(0)) 

'^ pr nt 

+ 2^ / / ^,{t- s)a^{u{s),u{t))dsdt^O. 



Consequently, recalling p.2p . v{r) — 0, u(0) = 0, < 7 = max{7i,72} < 1, and 
the fact that Ai are positive definite, we have 

p\\uir)r + {l-j)\\vml.<0, 

that implies u = a.e., and this completes the proof. D 

Remark 3. Theorem 2 also holds for the simplified problem (|2.1fp . see Remark 1 
and Remark 2. That is, with the assumptions in Theorem 2, there exists a unique 
weak solution for the simplified problem. 

4. Regularity 

Here we study the regularity of the unique weak solution of (|2.6I) . that is, a 
solution of (|2.8p - (|2.10l) . We explain the limitations for higher regularity in Remark 
4. We also prove higher regularity of any order of the solution of models with 
smooth kernels in Theorem 4. 

Corollary 1. If u° eV,v" eH,ge WliiO,T); Hr„), and f e L2{{0,T); H), then 
for the unique solution u of (|2.8p - (|2.10p we have 

(4.1) uGLoo((0,r);l/), 7i eLoo((0, T);iJ), ii G L2((0, T); y*). 

Moreover we have the estimate 

ll"llL=o((0,T);y) + \\u\\l^{(O.T)-H) + ll"llL2((0,T);y*) 
^^■^^ < C{\\u''\\v + \\V^\\ + MwUi^^TYMr,) + II/IIl.((0,T);H)}. 

Proof. It is known that if u„i -^ u, then 

|ii|| < lim inf ||wm||. 

Then, by p.l4p and the a priori estimates p.7p . we conclude (|4.ip and (|4.2p . D 

We note that, using Remark 3, Corollary 1 applies also to the simplified problem 

It is known from the theory of the elliptic operators, that global higher spatial 
regularity can not be obtained with mixed boundary conditions. Therefore we 
specialize to the homogeneous Dirichlet boundary condition, that is Fn = 0, and 
assume that the polygonal domain Q, is convex. We recall the usual Sobolev spaces 
jjr _ ffr^^-^ g^j^j ^g note that here V = -ffo(i7). We then use the extension 
of the operator A to an abstract operator with 2?(A) = H^{Vt)'^ n V such that 
a(w,w) = {Au^v) for sufficiently smooth u,v. We note that, the elliptic regularity 
holds, that is, 

(4.3) ||m||//2 <C||Am||, u^H^{VlYf\V. 
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Theorem 3. We assume that Tn — 0, and 
2 „t 



(4.4) 



> / Piis) ds < 1 or / max/3i(s) ds < -. 
t^Jo Jo '=i'2 2 



// m" G -ff'^, v'^ G y, anrf / G L2{{0,T); H), then for the unique solution u of 
(|2:8l) - (|2l0|) we have 



uGLoo((0,r);i/), mgX2((0,T);T/*). 
Moreover we have the estimate 

ll'"lli=o((0,T);H2) + ||u||L^((0,T);y) + I|m||l^((0,T);_H-) + 11 '^*' II L2((0,T);y*) 

<C{||w"|U. + ||^;"||y + ||/||Hl((o,T);ff)}. 
Proof. Differentiating p.Sp with respect to time, with notation v = v, we have 

P(iim(i)>¥'fe) + a(u„(t),<^fc) - V / ^.i{t-s)a^{u^{s),ipk)ds 
(4.7) , '-' 

= (/W,¥'fe) + £/3»Wa«("m(0),¥>fc), fc = l, ...,TO, tG (0,r), 



with the initial conditions 



Wm(0) = Mm(0) = "^{V^ , ipj)ipj , 

(4-8) T; 

^„(0) = ^„(0) = ^ (/(O) - Aur^iO),^,)^,. 
i=i 

Then, using j3i{t — s) ~ Ds^i{t — s) from p.ip and partial integration in time, we 
have 

2 

1=1 
2 „t 

+ y2 ^ii't^s)a,{u.^is),(pk)ds 
»=i-^o 

2 

= im,^k)+Y.(3^it)a^iu^{0),ipk) 
i=l 
2 

-5I^^W«(^™W,(^fe), fc=l, ...,m, iG(0,r). 

i=l 
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Now, multiplying by dk{t), summing k — 1, . . . , m, and integration with respect to 
t, we have 

p||i*™(i)f + (l-7)IKWIIy 
i^i Jq Ja 

< P\\u,niO)f + {I - l)\hnM\\v 

+ 2 / (/(r),^^(r))cir + 2^ / A(r)a,(u„(0), w^(r)) dr 

2 „t 

-2^] / 6Wa,(li,„(0),u„(r)) dr, 

where we recall that ^i(O) = 7i, 7 = max{7i,72}, and 7 = min{7i,72}. Then, 
recalling the fact that £,i are positive definite p.2p and integration by parts in the 
last term, we obtain 

<p\\Um{0W + {l-2)\\Uraiml 

+ 2 f{lir),u^{r))dr + 2Y, I P^{r)a^{u^,{Q) ^^{r)) dr 



^_^ 

2 „t 
-2^ / |i^{r)a,{u^{Q),u^{r)) dr 

2 2 

-2^^,(i)a,(u„(0),u„(t)) + 2^e,(0)a,(«„(0),u,„(0)), 

4=1 j = l 

that, using the Cauchy-Schwarz inequality, ||/3i||Li(E+) — 7i, 6(i) ^ ^i(O) = 7ii &n.d 
(|2.7p . implies 

p\\umit)r + {i-i)\hmmv 

+ 2/C1 max Ilit^Wf + Ci( / ll/WIMr)' 

0<r<t ^Jo ~ ' 

2 . . 2 ^ 

2 



2/^2(^7.) Ikm(0)||^.+2(^7.)C2 max ||ii„( 

i=l j=l ~ ~ 

2 2 

2/^3(5^7.) IK.(0)fy + 2(^7.)a3 max \\u^{r)fy 

2 2 2 

2M(E7.)ll^™(0)llv + 2(E^«)^4||u„Wfy + 2(5^7.) llli™(0)|| 



This implies, for some constant C — C(7i,72,p, T), 

lli^mllL^((0,T);H) + ll2^mllLoo((0,T);y) 

< C{||i*„(0)f + 1|U„(0)||2, + \\u^m?H^ + ll/llL((0,T);i/)}- 
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Then recalling u — u, the initial data from (j4.8p . and using 

||u™(0)||^. <||u°l|^., \\uUmv<\\v'>\\v, 
we have 

ll"™llioo((0,T);_H) + II"™IIloo((0,T);V) 

<C{\\uYh^ + \\vYv + \\fm? + \\l\\U(,,Ty.H))- 

We now find a bound for |jum(t)||jj2. We recall the eigenvalue problem (I3.3P 
with eigenpairs {{^iT'-Pi)}'fLi- Then we multiply (|3.5I) by Xkdk{t) and add for 
fc = 1, . . . , m to obtain 

2 ^t 



This, using the Cauchy-Schwarz inequality and (j2.7p . implies 
(4.11) 



\Au„M\' < l{\\f{t)f + p'\\u„M\^)+e\\Au„M\^ 



V / /3,(s) ds) max ||Au„.(s)f , 



that, by elliptic regularity (|4.3p and assumption (|4.4I) . gives us 

II"™IIloo((0,T);H2) < C'(^II/IIlo„((0,T);H) + ll^"lllL^((0,T);ff)j- 

From this and (14.91) we conclude 



ll"™llioo((0,T);H) + ll'^m|lLoe((0,T);V) + 11""™ llLo„((0,T);if2) 

< C{\\uYh^ + lk°fy + Il/IIL((0,T);//) + Il/IIL((0,T);«)}, 

that using ||/||L„((o,T);ff) < C'll/lli4'ii((o,T);ff)> by Sobolev inequality, we have 

ll'"™llioo((0,T);ff) + ll^ni|lL^((0,T);y) + 11^™ II L^ ((0,T);H2) 

<C{||u"|!l,. + ||r;"||Ull/ll?.H(o,T);/.)}- 
Finally from (14. 7p , similar to the proof of Theorem 1, we obtain 

\\'^m\\L2{{0,T);V) ^^{\W ll_H"2 + ||w llv + ||/||_H-i((o,T);H)|- 



The last two estimates then, in the limit, imply (|4.5p and the desired estimate (|4.6p . 
The proof is now complete. D 



Remark 4. If we continue differentiating (|4.7p in time to investigate more regu- 
larity, we obtain 

2 „t 

Pi'UmW^Vk) + a{Bmit),Vk) -Y^ / A(i-s)a-*(wm(s),'^fc)ds 



2 



= (/(i),¥'fe) + 5I/^»W«^("™(o)'^fe) 

2 
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Further, from $i{t)ai{u,n{0),fk), i = 1,2, we get $i{t)ai{um{0),u^{t)), but the 
Pi are not integrable. Besides, after integration in time, we can not use partial 
integration to transfer one time derivative from /3 to iim(^)j since /? is singular at 
t — 0. This means that we can not get more regularity with weakly singular kernels 
Pi. This also indicates that with smoother kernel we can get higher regularity in case 
of homogeneous Dirichlet boundary condition under the appropriate assumption on 
the data, that is, more regularity and compatibility conditions. 



Remark 5. For the simplified problem (|2.11l) . the inequality (14.111) is 

2 



\Au^{tW < -[\\fit)f + p'\\urn{tW)+e\\AuUm' 
(3{s) ds) max ||yl-u„i(s)||^. 







0<s<t 



Hence, the assumption (|4.4|) can be ignored, since /q/3(s) ds < j < 1. That is. 
Theorem 3 applies also to the simplified problem (|2.1ip . ignoring the assumption 

dUD. 



Remark 6. We recall the definition of the operators A, A\ and A2 from (|2.5p . and 
the fact that they are self-adjoint, positive definite linear operators. If the solution 
u is regular such that its second order partial derivatives are comutative, then 
the operator A^l"^ is comutative with the operators A-( , A^ ■ Here, the operator 
A^ {I e R) is defined by, see e.g., 03] 

00 

fc=i 

where {(A^, 'Pk)}kLi *re the eigenpairs of the operator A, and in a similar way A\ 
and A2 are defined. In this case the assumption (|4.4|) is replaced by 

t 
max/3i(s) ds < 1, 



since in (|4.10p we have 

2 pt pt 

y^ / Pi{t-s)ai{Uyn{,s),AUra{t))ds< 



max I3i{t - s)a{u,n{s) T Aum.{t)) ds 



t 

2 



< ( / max/3i(i — s) ds] max ||Aum(s)| 

\ /n i— 1,2 / 0<s<i 



where we used the fact that 

a.{v,Av)^{A^/^A]'^A'J\,A^/\) 

= {A]'^A^'^A\'\,A^'^v) = a,{A^/\,A^/^v) > 0. 

In the next theorem we prove regularity of any order of the solution of models 
with smooth kernels. 

Theorem 4. We assume that Tn = <2S, and condition (14. 4p holds. Assume (r — 
0,1,...) 

u"eH^+\ «"ei7^ ^eL2((0,T);iJ'-'=)(fc = 0,...,r), p,eWr^{0,T), 
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and the r -order compatibility conditions 



,0 „,o 



„,o 



u" e V, u'l := V 

«^=-(/(0)-^^i°)e^, 

P 



if r = 2 



r-3 2 



1 / d^~^ ' '^ ^ d^ 

P j=0 1=1 

Then for the unique solution u of (|2.8p - (|2.10l) we have 

■^u e L^((0, T); if '■+!-'=) (fc = 0, . . . , r + 1), 



and we have the estimate 

Ed M 



fc=0 



Loo((0,T);ff'-+i-'') 



\ ^^ II dt^ 

k=0 



L2({a,T);H^->') 



+ ||w°||^.+. + ||T;"||ff, 



We note that, with (3, e P^[~^(0, T) we have /3, G C^^^^ T] by Sobolev mequal- 
ity. Therefore u^ in the compatibihty conditions is well-defined. 

Proof. The proof is by an induction, the case r = following from Corollary 1 
above. 

Assume next the theorem valid for some ?Ti e N, and suppose 

7.0ei/''+2, i;OGi/'^+\ ^eL2((0,T);i/''+i-'^) (fc = 0,...,r + l). 

Suppose also the (r + l)*'^-order compatibility conditions hold. Differentiating the 
model problem (|2.6p with respect to t, we chack that u := -it is the unique weak 
solution of 

pu{t) + Au{t) 

2 .t 2 

- E / P^^^- •^)^«^(«) d^ = Lit) + E |3^{t)A^u" in fi X (0, T), 



u = Q 

li(-,0) = u°, u(-,0)=w° 



onrDx(0,r), 
in i7, 



for 



7.0 = t;", 



1 



i;"-ii(0) = -(/(-,0)-Aw°), / = /. 
P 

In particular, for r — we rely upon Theorem 3 to be sure that u G L2{(0, T); V), 
u' G L2((0,T);i/), u" G L2((0, T); y*). 

Since /, M°,t'° satisfy the (r + l)"^-order compatibility conditions, it is easy to 
see that f,vP,v° satisfy the r"^-order compatibility conditions. Thus applying the 
induction assumption for u we have 

d'^u 
IF 
with the estimate 



GL.o((0,T);iJ'^+i-'=) (fc = 0,...,r + l). 



''■+1 ,, jfe 
— ' d u 



dt^ 



fc=0 



L^CCCTji/f+i-") 



\ ^^ II difc 

fc=0 



L2((0.T);H'-'=) 



+ \\A\h^^^ + WAW 



18 



F. SAEDPANAH 



Since u = u we can write 




Y^ d^u 

^ df" 
fc=i 


Loo((0,T);_H'-+2-fc) 




L2((0,T);H'-+i-'=) 




L2((0,T);H'-+i-'') 



||«°||ff,.+i + P«°|lH. + ||/(0)||ff,. 



l"°llff-+= + ll«°lk.-+l 



•fc=0 



where, for the last inequahty, we used the Sobolev inequahty 

|l/llc([0,T];ff'-) < C'(||/||l2((o,T);H'-) + ||/||L2((0,T);ff'-))- 



We now find a bound for 



a{u,A''+'^u)^{f~pu. 



|w||L^((0,T);if+ 



2 . We set w = A'' 



in (US to obtain 



2 pt 

A'-+iu) + V / (5,(t ~ s)a,{u{s),A'-+^u{t))ds. 

,-1 "'0 



This, using the Cauchy-Schwarz inequahty and p.7p . imphes 



(tl>' 



s) ds] max IIA 2 M(s)|h 

' 0<s<i 



that gives us, by elhptic regularity (|4.3p . 

II"IIlo„((0,T);H'-+2) < C'(^II/IIl^((o,t);H'-) + II^IIloo((o,t);H'-; 
Adding this inequality to (|4.12p we deduce 

r+2 



E\id''u 
lldi^ 
k=0 



L^((0,T);ff'-+2-fc) 



r+1 



\ ^^ II dt*^ 

fc=0 



L2((0,T);_ff '■+!-'=) 



|u"||^,-+. + ||z;0||h.+i 



that is the assertion of the theorem for r + 1. Now the proof is complete. 



D 



We note that Remark 5 also holds for Theorem 4. Remark 6 can be applied 
to Theorem 4, provided the solution u is smooth enough such that the operator 
A'^ is comutative with the operators A^' , i — 1,2, that is, when (r + 2)-th order 
partial derivatives of the solution u are comutative. 
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